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We propose a new mechanism for baryogenesis at the 1-200 MeV scale. Enhancement of CP 
violation takes place via interference between oscillations and decays of mesinos —bound states of 
a scalar quark and antiquark and their CP conjugates. We present the mechanism in a simplified 
model with four new fundamental particles, with masses between 300 GeV and 10 TeV, and show 
that some of the experimentally allowed parameter space can give the observed baryon-to-entropy 
ratio. 


I. INTRODUCTION 

Baryogenesis, the explanation for the asymmetry be¬ 
tween matter and antimatter in our universe, is a pro¬ 
found puzzle for particle physics, as well as one of the 
strongest motivations for extending the standard model 
(SM). The first proposed solution, by Sakharov in 1967, 
laid out three conditions for a successful resolution to the 
puzzle [T . The first condition, C and CP violation, is sat¬ 
isfied by a CP-violating phase in the CKM matrix of the 
SM, however the effects of this phase are too suppressed 
in the early universe to produce a sufficient asymmetry. 
The second condition, violation of baryon number, is sat¬ 
isfied by anomalous electroweak processes which are suffi¬ 
ciently fast at high temperature to produce baryon num¬ 
ber [21, but only if the third condition, departure from 
thermal equilibrium, is satisfied. The minimal SM with a 
Higgs boson at 125 GeV does not have any phase transi¬ 
tion or sufficiently long-lived heavy particles to produce a 
sufficient departure from thermal equilibrium for baryo¬ 
genesis PJ. 

In most supersymmetric extensions of the standard 
model, cosmology theory favors a low reheat scale af¬ 
ter inflation Furthermore, in many theories, the 

absence of observed isocurvature perturbations favors a 
low inflation scale Eli. A low inflation or reheat scale 
would imply that the baryogenesis scale must also be 
low, in some cases as low as an MeV. Although propos¬ 
als for baryogenesis at scales as low as an MeV exist [7], 
the time scale at this epoch is relatively long compared 
to typical particle physics scales, making departure from 
thermal equilibrium difficult to arrange. Models which 
do depart from thermal equilibrium at these low scales 
then produce additional entropy which dilutes the baryon 
abundance. In addition, CP violation at these low scales, 
if large enough to produce enough baryon number to sur¬ 
vive the dilution, can easily lead to particle electric dipole 
moments which are in conflict with experimental bounds. 
Finally, low energy baryon number violation can lead to 
rare decays in conflict with observation. It is therefore 
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difficult to produce enough baryons at low energy while 
satisfying experimental constraints, motivating the ex¬ 
ploration of new mechanisms for low scale baryogenesis, 
particularly ones which enhance the quantum mechanical 
interference which is necessary for CP violation. 

In the SM the CP-violating phase is unphysical if any 
of the small mixing angles or like-sign quark mass differ¬ 
ences vanish. A reparameterization-invariant measure of 
CP violation |FJ in the SM is thus very tiny. Neverthe¬ 
less, the SM does provide a wealth of large CP-violating 
asymmetries in the decays of oscillating neutral mesons, 
illustrating the efficiency of oscillations for enhancing the 
effects of CP-violating phases. A proposal [3] that a sim¬ 
ilar enhancement of CP violation could occur during an 
electroweak phase transition was shown not to work due 
to the too rapid rate of thermal decoherence [TO]. 

In this paper we present a new baryogenesis mech¬ 
anism. We show how in an extension of the SM, 
the particle-antiparticle oscillations of mesinos mm 
(bound states of a fermion quark and a scalar antiquark 
or vice versa), when combined with baryon-violating 
scalar decays, can give rise to baryogenesis. The scalar 
quark, whose mass is of order a TeV, is produced by the 
out of equilibrium decays of a long lived heavy neutral 
fermion, below the QCD hadronization scale but before 
nucleosynthesis. Both the oscillation rate and the de¬ 
cay rate of the mesinos are rapid when compared with 
the decoherence time, allowing quantum mechanical in¬ 
terference between particle and anti particle decays to 
produce substantial CP violation. This mechanism for 
producing baryons is amusingly similar to the mecha¬ 
nism for CP-violating production of charged leptons in 
the neutral kaon system. We will illustrate the mecha¬ 
nism in a simplified model which introduces a minimal 
number of new particles beyond the SM. Our mechanism 
and the new particles could be part of a supersymmet¬ 
ric extension of the SM with baryon-number-violating 
R-parity violation. 

The organization of the paper is as follows. In Sec. [IT] 
we introduce the model, describing the CP-violating os¬ 
cillations of the mesino. Section |III| details experimental 
constraints on the model, from production at colliders 
and precision measurements of rare processes. In Sec. |IV| 
we describe the cosmology of the model, estimating the 
parameters needed to explain the observed baryon asym¬ 
metry. We conclude in Sec. [V] 
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II. THE MODEL AND SOME PHENOMENA tive parts, just like the case of the neutral mesons, 


The ingredients we need for our model are a complex 
scalar </> and n Weyl fermions JV,. (f> is a color triplet, 
SU(2)i singlet, and carries hypercharge —1/3. Each IVj 
is taken to be a singlet under the SM gauge group. The 
relevant interactions of these are 

H, _.) y i: j(f>diNj —TTijyijNiNj T otjj& d t Uj T c.c. (1) 

Ui and di are the left handed up and down-type singlet 
antiquarks, respectively, with the subscript labeling the 
generation. Another possible baryon-number-violating 
interaction of (j) , (p q t q r , where the q are quark doublets, 
is neglected because we are considering only the inter¬ 
actions that could come from baryon-number-violating 
interactions in a supersymmetric theory where <fi is a su¬ 
perpartner of a down-type singlet quark. By rotating and 
rephasing the IVj we can make the singlet mass matrix 
toat real and diagonal, atj is a 3 x 3 matrix containing 
nine complex parameters, with seven phases removable 
by reparameterizations. The remaining phases in this 
matrix do not play a role in our baryogenesis mechanism 
and will be ignored for the rest of the paper. Having 
exhausted our freedom to rephase the fields, each of the 
elements in the 3 x n matrix contains a physical, CP- 
violating imaginary component. 

We will see that the simplest version of the model in¬ 
volves three new Weyl fermions, i.e. n = 3, which we 
label in order of their masses, itin 1 < itin 2 < toat 3 - The 
colored scalar is produced in the early Universe by late, 
out-of-equilibrium decays of IV 3 . The scalar forms color- 
singlet bound states with SM quarks, termed mesinos, 
which undergo particle-antiparticle oscillations and de¬ 
cay. The two lighter singlets provide common intermedi¬ 
ate on- and off-shell states for the mesino and antimesino 
that allow for them to violate CP (and baryon number) 
in the interference between decays with and without mix¬ 
ing, sourcing the baryon asymmetry of the Universe. 

In the following sections we discuss mesino-antimesino 
oscillation in this model and the constraints from collider 
experiments and precision measurements. 


A. Mesino oscillations and decay 


h = m- 1 - T. (2) 

The eigenstates of H, with eigenvalues 0Jl,h, can be 
given in terms of the flavor eigenstates |<f> g ) and |$ g ), 

|*L.ff>=p|^>±#,>, (3) 

with L and H referring to light and heavy. The complex 
numbers p and q are related by 

fo\ 2 _ Mf 2 -(»/2)n 2 
W m 12 - (i/ 2 ) r 12 ■ 

The mass and width differences between |$a) and 1$#) 
are 

Ato = ttih — tul = Re (w# — ujl) , (5) 

AT = T h -T l = -2Im(wi/-Wi), (6) 

with 



In the scenario we will study, the mass difference is small 
compared to the masses of the heavy and light eigen¬ 
states, 


m H + rn L 

m H — m. L ^ ---= . 


( 8 ) 


Furthermore, the mass of the mesino is mostly supplied 
by the (j >, ~ m^. As in the case of meson oscillations, 

it is often useful to introduce a dimensionless parameter 
to characterize the oscillation rate, 


A TO 


x = 


(9) 


where T = (r//+Pi)/ 2 is the average lifetime of the 
states, x 1 indicates that the mesino system oscil¬ 
lates rapidly before decaying while x <C 1 means that 
the mesino typically decays before oscillating much. 

We use | ( t )) to label a state at time t that began at 
t = 0 as a pure mesino and \$ q ( t )) to label one that be¬ 
gan as a pure antimesino. These evolve in time according 
to 


At temperatures below the QCD hadronization scale, 
T c ~ 200 MeV, if the colored scalars are sufficiently long- 
lived, they will bind with light SM quarks to form mesinos 
and antimesinos. In analogy with the naming conven¬ 
tion for mesons, we will refer to the (electrically neutral) 
mesino containing </>* and q = d, s,b as with its an¬ 
tiparticle containing <f> and q = d,s 1 b. 

and d> g form a pseudo-Dirac fermion and, as dis¬ 
cussed in Ref. m, the system can be described using a 
two-state Hamiltonian containing dispersive and absorp- 


l$g (t)) = 9+ C t ) l $ g> - Jj9- (t) ( 10 ) 

\$ q (t))=g + (t)\$ q )-^g-(t)\$ q ), (11) 

with 

g± (t) = i ( ' e -irn H t-±T H t ± e -im L t-\r L t^ _ ( 12 ) 

The baryon-number-violating operator (ff dplj and its 
conjugate in Eq. ([l]) allow for the mesino and antimesino 
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flavor eigenstates to decay to collections of hadrons with 
baryon number B = +1 and B = —1, respectively, with 
amplitudes related by complex conjugation. In terms of 
the time-independent amplitude for a mesino to decay to 
a state with B = +1, A4 , the time-dependent amplitudes 
for initial mesino or antimesino states to decay to B = ±1 
(denoted B and B ) are 

{B\$ q (t)) = g+ ( t)M , 

(B\* q ( t)) = - q -g- V)M*, 

p ( 13 ) 

{B\* q {t)) = -~ q 9- (t)M , 

(B\* q ( t))=g+ 

Squaring these and integrating over t, we can find the 
probability that an initial mesino state decays to B and 

B, 


pOO 

P^BOC / dt\(B\$ q (t))\ 2 , (14) 

Jo 

p oo 

P^Btx dt\(B\<S> q (t))\\ (15) 

Jo 

and similarly for an initial antimesino. Using these ex¬ 
pressions, we write down the baryon asymmetry per 
mesino-antimesino pair, 


e_e = Ab x Br^^s, 

where 

. ^ + P$> q ^B ~ P$ q ^>B 

A.b = -• 

+ P^ q -yB + P® q ^B + P$ q -yB 


(16) 


(17) 


and Br^^B = r$ 9 _s.B/T is the (time-independent) 
branching ratio of a mesino flavor eigenstate to decay 
into a state with B = +1. In terms of the elements of 
H , Ab is 


2 ImIW 1 * 2 r 12 
r 2 + 4|M 12 | 2 ' 


(18) 


To get a sense of how large the asymmetry es can be, 
we estimate the magnitudes of the elements that appear 
in it below. 


1. Estimating M 12 


The leading contributions to M 12 in the system 

arise from the diagrams in Fig. [l] resulting in mm 


M 12 = ^M 12 (Ni) 

i 



^ m\r — m 

7 


2 

<E> 


q 


(19) 


where we have written the contribution from each Ni as 

Af 12 (Ni). 



FIG. 1. Tree level contributions to M\ 2 - 



The mesino decay constant /$ can be related to the 
B meson decay constant EL 


/$, — /b. 



a a (mb) 

da (m t ) 


6/23 


/ a s (mt) 
Vo; s (m^) 


6/21 

■ ( 20 ) 


The cosmology of the model will motivate focusing on the 
strange mesino; to estimate its decay constant, we use 
/b s = 225 McV [I5j, MS quark masses, mb = 4.18 GeV, 
mt = 160 GeV jTB], and leading order QCD running (the 
contribution of which is negligible above the top mass in 
this case), and find 


U s ~ 21.5 MeV< 


1 650 GeV 

TO 0 


( 21 ) 


In what follows, we will take q = s, specifying to the case 
of the strange mesino, 4v 


2. Estimating F and Ti 2 


As mentioned above, the operator <fi*diUj in Eq. ([lj| 
allows for the mesino to decay to a collection of hadrons 
with B = +1 with a rate 

r^-s-s = \ a n\ 2 ( 22 ) 

i,3 


where we assume that we can ignore the masses of the 
hadrons in the final state in comparison to the mesino’s 
and we have defined a 2 B = JV • | | 2 . 

In general, the mesinos can also decay to a singlet plus 
hadrons, 4> g —»• Nj +h where h labels the hadronic state, 
through the operator 4>diNj with a rate proportional to 


y.jj | . This is kinematically allowed for m$ q > my +m^. 
If h is self-conjugate under C (h = tt°, , p,...) then 
this is a final state for antimesino decay as well, leading to 
a nonzero value of Ti 2 and CP violation through Eq. (181. 
This CP asymmetry is suppressed when |Fi 2 | /T is small; 
for this reason we take only TV! to be kinematically al¬ 
lowed in the decay of the strange mesino, and set its 
mass close enough to the mesino’s so that single hadron 
final states dominate the partial width. In this case 
\r 12 \ ~ T^jVi an d |ri 2 |/r can potentially be unsup¬ 
pressed. 

Assuming that ~ mjv 4 , we can approximate T 12 
as being dominated by the decays <4> s , <E> S —► Ni+rj, where 
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FIG. 2. 2-body contribution to Ti 2 from <F S , < 1> S —>■ IVi + r? 
where r/ is the lightest meson containing ss. 


ij is the lightest pseudoscalar meson composed of ss, as 
shown diagrammatically in Fig. [2] This results in 


r t-i2 —body 

12 — 1 12 

yh^N, 


167T 


|F(u4)| 2 A 1/2 



(23) 

(24) 


where A (a, b, c) = a 2 + b 2 + c 2 — 2ab—2ac—2bc and F ( q 2 ) 
is the 4> s —> P transition form factor at momentum trans¬ 
fer q. We provide estimates of the form factor F (m ‘ 2 Nl ) 
in Appendix [A] For m$ s ~ 1 TeV with a splitting 
~ 1 GeV, the form factor is O (lO -3 — 10 -2 ). 
Defining the mass splitting, A= |nT$ a — and 
assuming that they are small compared to the mesino 
mass, the rate to decay to the lightest singlet and the 
off-diagonal term in the absorptive part of the Hamilto¬ 
nian are given by 




rJVi 


~ |r 12 | ~4x 10" s GeV 


y s i 


F (m 2 Ni ) 


10 


-2 


0.1 

( Am, 


4>N i 


V 1 GeV 


(25) 


choosing to normalize the mass of the scalar above exper¬ 
imental constraints from collider searches (see Sec. Ill A 
below). Comparing this with Eq. (251, we expect 


that, barring extreme cancellations between contribu¬ 
tions from the different singlets, |Mi 2 | Ti 2 , 

The asymmetry per mesino pair can be written sim¬ 
ply in terms of the oscillation parameter x in this limit. 
When \T 12 \ <C |Mi 2 |, x ~ 2 |M i2 | /T and 


x r sin /? |ri 2 | 

£B = ~T“i— 2~~ 

1 + x z I 9 


■O 


M 


12 


(27) 


where 


M 12 (Vi) 

Mi 2 


0 < r < 1 , 


(28) 


and the reparameterization-invariant, CP-violating 
phase is 


0 = arg{[M 12 -M 12 (IVi)]*r 12 }. (29) 

/3 is defined this way so as to be generically O { 1 ), making 
use of the fact that Ti 2 and (W) have the same 
phase. If we also make use of the relation between the 
decay rate to a singlet and the off-diagonal term in T, 
|r i2 | ~ r 4 , g ^ J Vi = TBrj^jVj then 

xr sin/3 

es — ——;— 5 ~Br$ _>.BBr$ (30) 

1 + x z 

The question of whether 2 is large or small com¬ 
pared to the total mesino decay rate T is determined 
by the partial width to antibaryons, T$ If a B is 
much smaller than about 10 ~ 4 |y sl | (650 GeV/m^), then 
\M 12 \ » and i»l. Then, 


for a scalar of mass ~ 1 TeV. 

For the mesino to hadronize before decaying, the de¬ 
cay rate has to be slow compared to the timescale for 
hadronization, F -C 1/fhad ~ Aqcd ~ 200 MeV. The 
partial width to satisfies this condition for any set 
of parameters we consider. Requiring that the partial 
width to baryons is smaller than Aqcd limits a B < 
0.1^650 GeV/m^. 


e B -sin/3Br$ _>.BBr$ -yjvj. 

X 


(31) 


In this case, by inspecting Eqs. (26) and (251, typically 


-kb- 


x < 10 2 , depending on the hierarchy of |ri 2 ]and 
The asymmetry is suppressed by x and a potentially 
small branching ratio. On the other hand, if a B 
10 ~ 4 |y s i| (650 GeV/m 0 ) then |T 12 | < |M 12 | < 
and x -C 1. The asymmetry becomes 


3. Simplifying the asymmetry 


e B -X x r sin /3 Br^^j^, (32) 


In the limit of small mass splittings between and 
IVj, ATO 0 JV; then the contribution to Af 12 from 

Ni is approximately 

|M 12 (iVi)| ~ 1.54 x 10" 6 GeV|^-| 2 

/ 650 GeV A / 1 GeV \ ^ 

V m cj> ) V ^ m 4>Ni ) ’ 


suppressed by the small values of both x and Br.j, ^^,. 

The asymmetry is typically largest when T<i, 9 ^. b ~ 
| r 12 1 < |Af 12 |. In this situation Eq. (31) applies and 
e B ~ 10 _ 2 rsin/3. We note that it is possible to achieve 
an asymmetry as large ases = l/ 8 if T.j^b — |Ti 2 | — 
M-\ 2 1. However, this requires large, unnatural cancella¬ 
tions (at the level of 10 ~ 2 ) between the contributions to 
Mi 2 . 
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III. EXPERIMENTAL CONSTRAINTS 
A. Collider Constraints 


In this model the colored scalar can be pair produced 
in pp and pp collisions which then decay to multijet fi¬ 
nal states. If the dominant decay mode is through the 
coupling aij(f)*diUj + h.c., then <j> will decay to two hard 
jets, leading to a four-jet signature. If the decay of <f> 
proceeds mostly through Ni, the situation is somewhat 
more nuanced as we describe below. 

In proton-(anti)proton collisions, the cjxj)* pair is pro¬ 
duced in the lab frame with very little boost in the di¬ 
rection transverse to the beam. The <j> decays to a very 
soft light quark with a momentum of O and iVi 

that is also essentially at rest in the (f> rest frame, because 
of the small mass splitting between the (j> and Ni. The 
N\ decays through an off-shell </> to three quarks with a 
rate 

dT Nl \yn\ 2 a B P ( m Nl ~2p\ 2 

dp “ 5127T 3 ttin 1 Vjj + ) 


ignoring the masses of the quarks, p is the momentum 
of the quark produced in the IV! —> (f>di splitting in the 
IVi rest frame. Because of the small mass splitting, the 
momentum of this quark is peaked at small values—in 
most decays p < yj A niNy ■ This means that, for a 
splitting A = 1 GeV, N\ masses below ~ 1 TeV 
result in one of the three quarks in IVi decays having a 
transverse momentum in the lab frame of pt < 30 GeV, 
making them hard to observe. Therefore, at colliders, </>’s 
either decay directly to two hard jets or appear as two 
hard jets in their decays to iVi. 

To use collider searches to set limits, we will fo¬ 
cus on the situation where the scalars are produced 
and decay promptly. At the Large Hadron Collider 
(LHC) this requires decay lengths cT ( f,,CTj^ 1 < 1 mm. 
If the </>’s decay directly to two quarks this requires 
as > 10”V 650 GeV/m^. If, instead, (j) decays pro¬ 
ceed through Ni, then (£V \yn\ 2 Y^ 2 > 10 ” 4 , for a 
form factor of 10 ” 2 and a mass splitting of 1 GeV 
and | m| 2 ) 1 / 2 a fl > 10 _6 a/650 GeV/mjVy■ We will 
see in Sec. IIIB that \yn\ is essentially unconstrained 
so this last constraint can be taken to be as > 
10 _6 \/650 GeV/mjVj. The collider limits we will outline 
are relatively insensitive to the particular values of these 
couplings so long as they are large enough for prompt 
decays, given <p(f>* production is dominated by QCD. We 
can also safely assume that these couplings are not so 
large as to make the width of (/> large enough to impact 
the limits. 

The requirement of prompt decays at hadron colliders 
is not a strict condition for the viability of the model, but 
merely simplifies our analysis later; a more robust lower 
bound, weaker on the couplings by O ( 10 6 ), comes from 
not spoiling Big Bang nucleosynthesis (BBN). Indeed, 
long-lived particles at colliders have been studied as a 


probe of baryogenesis, see, e.g., Ref. da. In the long- 
lived case, the lower limits on the particle masses can be 
quite strong, as high as ~ 1 TeV m- 

A search for squarks decaying to a b and an s quark 
was undertaken at the LHC by the ATLAS collaboration 
using 17.4 fb ” 1 of 8 TeV pp collisions |T5j. Squarks with 
masses between 100 GeV and 310 GeV were ruled out at 
the 95% confidence level (CL). Since they are both color 
triplet scalars, the production cross section for a <p<p* pair 
is equal to that of a squark-antisquark pair. Thus, this 
exclusion would directly apply to our model if the leading 
ctij were « 3 i or 0 : 32 , and constrains m$ > 310 GeV0 In 
addition, ATLAS has searched for the pair production of 
scalar gluons that each decay to two gluons using 4.6 fb ” 1 
collected at 7 TeV [20], finding no excess and setting an 
upper limit on the cross section for this process. This 
limit applies to </> in the case that it decays to a pair of 
light flavor quarks, which is the case if the a t j , i. j = 1 , 2 
couplings dominate. Using NLL-f ast m to calculate the 
cross section to produce a cjxj)* pair at next-to-leading 
order with next-to-leading logarithmic corrections, this 
search sets the lower limit > 275 GeV at 95% CL. 

The CMS collaboration performed a search with 
19.4 fb ” 1 of 8 TeV data, looking for the pair produc¬ 
tion of resonances that decay to two jets [22] • Col¬ 
ored scalars decaying to a 6 and a light flavor quark are 
ruled out at 95% CL for 200 GeV < < 385 GeV 

and those decaying to two light quarks are excluded for 
200 GeV < < 350 GeV at 95% CL. 

Although we expect the tfr to decay to essentially two 
hard jets, as explained above, if it decays through IVi, 
there could be apparent six-jet events from (jxj)* produc¬ 
tion and decay. Searches for three-jet resonances [23] 
generally have more reach than those for two-jet reso¬ 
nances because of the comparative ease of distinguishing 
this from backgrounds. In cases where these decays in¬ 
clude heavy flavor, limits as high as 600 GeV can be set. 

In light of this discussion, we consider 600 GeV as a 
conservative lower bound on the mass of deferring a 
complete study of the collider limits for future work. For 
this reason, when specifying m^, we have taken 650 GeV 
as a benchmark value. 

The prospects for future searches for (j) —> 2 j are 
promising, with a study of the 14 TeV LHC indicating 
that uifp up to 750 GeV could be probed with 300 fb ” 1 
and up to 1070 GeV with 1000 fb ” 1 [25 . Such searches 
would push significantly into the region of parameter 
space compatible with an explanation of the baryon 
asymmetry of the Universe, as we will see in Sec. m 


1 A similar limit ~ 300 GeV was derived in Ref. m for mesinos 
oscillating and decaying to a top and a light quark, leading 
to a same-sign dilepton signal which was constrained by CMS 
data m- 
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FIG. 3. Diagram responsible for the dimension-9 A B = 2 
operator in Eq. (341 that leads to neutron-antineutron oscil¬ 
lations. 


B. Constraints from Rare Processes 

Our model possesses new interactions which violate 
CP, quark flavor, and baryon number, potentially leading 
to observable conflicts with Standard Model predictions. 
In examining the consequences, we note that our model 
could be part of an R-parity-violating supersymmetric 
model |25fB7| . and we could use some of the existing 
studies of constraints on such models [23H3I]. The con¬ 
straints on baryon-number-violating parameters in su¬ 
persymmetric theories depend on many soft supersym¬ 
metry breaking parameters which play no role in baryo- 
genesis. In this section we consider the constraints on 
the parameters in our more minimal model. Although 
baryon number is violated, given that that lepton num¬ 
ber is conserved 0 and given that the only fermions lighter 
than the proton are leptons, the proton is stable. There 
are stringent constraints on couplings involving the first 
generation from neutron-antineutron oscillations, and on 
baryon-violating couplings of the strange quark, from 
heavy nuclei decays. 

Dimension-9 six quark A B = 2 operators arise in the 
low energy effective theory from the diagram shown in 
Fig-0 The relevant operators are of the form 

( v~' yikUjk \ cq m a no - - - _ - _ . . 

> yry- - j — didjdiu m d n u 0 . (34) 

\k Mn * j < 

Constraints on such operators with various color and 
spin combinations from neutron-antineutron oscillations 
have been considered in many previous works HE ED HO 
I32H3U- An accurate computation of the rate requires 
perturbative renormalization group scaling and match¬ 
ing to lattice QCD matrix element computations of the 
various operators [33, [35], and currently still has uncer¬ 
tainties of order 1. Currently the best constraint on the 
neutron oscillation rate is 2.9 x 10” 33 GeV |34| . If we 
simply estimate the matrix element of 

(n\dddduu\n) ~ 0( 1(T 5 GeV 6 ) , (35) 


or conserved mod 2, if there are Majorana neutrino masses 


which is consistent with the lattice computations, we find 




-T<2.9x 10~ 28 GeV -5 

mr. 


(36) 


Neglecting the very weakly coupled iV 3 , and using ATy, ~ 
Mn 2 ~ to ^ ~ 650 GeV gives 

(l&+l&)“ii<°( 10 ~ 14 )- (37) 

The best constraints on operators involving strange 

quarks arise from limits on dinucleon decay into kaons jD 
HD OSH EMI ■ Using the rough estimates from Ref. [3U] for 
the nuclear matrix elements and the experimental bound 
from Ref. [36] gives a bound on the coefficient of the op¬ 
erator ssdduu which, for M R / 1 ~ Mn 2 ~ to^ ~ 650 GeV, 
translates into bounds on combinations of couplings 

(Usi + Vsi) a n < (38) 

(2/Si + Z/d 2 ) < O(10” 14 ), (39) 

(ydiVsi + VdiVsi) ai 2 a n < £>(10~ 14 ), (40) 


although the strong interaction and nuclear physics un¬ 
certainty in these bounds spans several orders of magni¬ 
tude. Neutron oscillations and nucleus decay give much 
weaker constraints on coefficients of operators containing 
heavier quarks, as the required loops involve weak inter¬ 
actions and small CKM parameters. See, e.g., Ref. m 
for a summary of such constraints in R-parity-violating 
supersymmetry. 

We also need to check that the new sources of flavor 
violation will not violate the stringent constraints com¬ 
ing from the agreement between theory and experiment 
in neutral meson oscillation phenomenology. After in¬ 
tegrating out 4 > and the Vj, the interactions in Eq. ([T]| 
will lead to A F = 2 flavor-changing operators.The most 
stringent constraints come from the neutral kaon oscilla¬ 
tions. 

The relevant AS = 2 Hamiltonian that is generated 
reads [3T 


n 


AS=2 

<f>Ni 


E 


y* di ydjysiy* sj 
2 9 3 3 7t 2 to| 


d a R rs a R d p R 7„ S 


FO 


(41) 


where a and f3 are color indices, and we have assumed 
that TYifj) ~ rriNi- This contributes to the K°-K° mass 
difference, 


A m K = 2Re(K Q \n^= 2 \K°) 

0 y dl ydjysiy SJ fukSk 
= Re S 283 4 T C m | ■ 

ij t 


(42) 

(43) 


with f R the kaon decay constant. Requiring that this not 
exceed the measured mass difference of (3.484 ± 0.006) x 
10~ 12 MeV [16] leads to a constraint on the couplings, 

/ \ i/4 


2 


( 44 ) 
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where we have used fx = 155 MeV [35]. Since the 
Yukawa couplings are complex, there is also a contribu¬ 
tion to CP-violation in the neutral kaon system. This is 
typically parameterized by 


lm(K°\'H AS= 2 \K 0 } 
\/2A rriK 


(45) 


If we assume that the kaon mass difference is saturated 
by the SM contribution, then the operator in Eq. (41) 
gives a contribution to CP-violation of 


/ 650 GeV \ * * , , 

e K = U[ - im 2^y di y d jysiy S j- (46) 

V m 4> J itj 

Requiring that this is less than the measured value, 
(2.228 ± 0 . 011 ) x 10~ 3 US], results in a constraint, 


the baryon-to-entropy ratio via t]b = 3.9 x 10 -9 x f Ibh 2 . 
Planck has measured fit,h 2 = 0.0221 ± 0.0003 [40, which 
gives tib = (8-6 ± 0 . 1 ) x 10 ” 11 . 

In the minimal version of the model, the colored scalars 
can be produced at late times, well after they have frozen 
out, through N :i decays. Those produced in decays af¬ 
ter the universe has cooled below the QCD hadroniza- 
tion temperature T c — 200 MeV bind immediately to 
form mesinos and antimesinos. The CP-violating oscil¬ 
lations and decays (to antibaryons and baryons) of the 
(anti)mesinos give rise to a baryon asymmetry. To quan¬ 
tify this asymmetry, we must first determine the number 
of N 3 present after the universe has cooled below T c . 

There are two processes which control the number den¬ 
sity of N 3 : annihilation to quarks and scalars and decay 
to a scalar and quark. The number density of N 3 at time 
t is then 


/ \ !/ 4 

pY *»«»•** j < “VeBSW (47) 

Since we expect that the CP-violating phases of the 
Yukawa couplings have no reason to be suppressed, the 
limit from ex is stronger than the one from A nix ■ How¬ 
ever, this constraint does not impact the viability of the 
model; even 0 ( 1 ) values of \y S i\ are allowed so long as 
12/di I < O (10” 2 ). 

Limits from B°-B° and B s -B s mixing are no more con¬ 
straining than from kaon mixing, and values of \ybi\ as 
large as unity are allowed. 

Experimental upper bounds on Electric Dipole Mo¬ 
ments (EDMs) place strong constraints on new CP- 
violating physics. For new CP-violating physics at the 
weak scale, there may be nonstandard contributions to 
EDMs at one or two loops, which place constraints on 
combinations of the new couplings and phases. The 
relevant constraints on R-parity-violating supersymme¬ 
try are summarized in Ref. m, and on supersymmetric 
models in general in Ref. |39| . In our minimal model, 
which does not include new couplings to leptons or to 
left handed quarks, and which does not require any new 
couplings to be large, there are no one loop contribu¬ 
tions or large two loop contributions to EDMs and no 
constraints on the CP-violating phases. 


IV. COSMOLOGY AND THE 
BARYON-TO-ENTROPY RATIO 


The baryon asymmetry of the Universe is quantified 
by the ratio of the net baryon number to entropy, 


Vb = 


n B 

s 


(48) 


The most precise determination of this quantity comes 
from measurements of the CMB that fix the baryon den¬ 
sity in units of the critical density, fib , which is related to 


—n relic -r N3 t 
UN 3 — n N 3 e 


ilrelic 

at 


(49) 


Here, n™' IC and a re i; c are the number density of N 3 and 
the scale factor, respectively, at T ss T ~ mjv 3 /20, and at 
is the scale factor at t. Similarly to N\ and IV 2 , we assume 
that N 3 is close in mass to q J, Ato^jv 3 -C TO 0 ,TOat 3 ■ Then 
its decay rate is given by 


r iv 3 - 

q—d,s,b 


y 2 q 3 Am lN 3 

8 ir m$ q 


(50) 


We want an appreciable number of N 3 to survive until 
T c ~ 200 MeV which corresponds to a time t c ~ 10 -5 s 
assuming standard thermal history. Our thermal history 
is not standard since we have N 3 decaying and inject¬ 
ing energy into the plasma competing with the cooling of 
the plasma due to the expansion of the universe. How¬ 
ever, as we show below, that does not drastically affect 
the time vs. temperature relationship. Therefore, to get 
a respectable number of N 3 to decay after T c , produc¬ 
ing mesinos, we must have Tjv 3 < l/t c ~ 10 -19 GeV. 
For a splitting A m^x 3 ~ O(GeV), this requires y 2 3 < 
10~ 15 (mAr 3 /TeV). Since the N 3 annihilation rate is pro¬ 
portional to y 4 3 , Yukawa couplings that are this small 
mean that the annihilation rate was always much smaller 
than the expansion rate of the Universe. Hence, assum¬ 
ing reheating happened at a high temperature and the 
N 3 were in equilibrium with the plasma, n 1 ^ 0 = (3/4) n 7 
where n 7 is evaluated at T = 771 x 3/^0 and only the decay 
rate F^ dictates the N 3 number density for T <T C . 

We can determine the baryon-to-entropy ratio by co¬ 
evolving the N 3 and radiation energy densities, 

df> ^ = -4 Hp lad + T N 3 m N 3 n N3 , (51) 

= -3Hp Ns - T N 3 m N 3 n N3 , (52) 

along with the baryon (minus antibaryon) number den- 











sity, 


= -3Hn B + ^AT N 3 e B n N3 . (53) 


Equations (511 and (521 describe the radiation and N 3 


energy densities, taking into account the expansion of the 
universe and the fact that the decays of N 3 heat up the 
plasma. Equation (53I describes the evolution of the net 


baryon number density, which develops a nonzero value 
due to CP-violating oscillations of mesinos produced by 
N 3 decays. The factor of 1/2 appears due to the def¬ 
inition of e B hr Eq. (16) as the baryon asymmetry per 


mesino-antimesino pair. A counts the fraction of mesino 
states that undergo CP-violating oscillations. In this 
work, we focus on the strange mesino since, being an 
isoscalar, it is not decohered by scattering on the pions 
present in the plasma after confinement, a possibility in 
the case of the down mesino that would complicate the 
analysis. Thus, we take A = 1/3, assuming that equal 
fractions of up, down, and strange mesinos are formed 
which should be true to ~ 30%, the level of SU(3) flavor 
breaking observed elsewhere. It is important to keep in 
mind that in additon to Eq. (531, the net baryon num¬ 


ber density is subject to the constraint that n B = 0 for 
T < T c since scalars produced above the QCD confine¬ 
ment temperature do not hadronize and therefore do not 
undergo coherent oscillations. 


Instead of solving Eqs. (511 to (531 exactly, we first 


use the sudden decay approximation to gain an intuitive 
understanding of the baryon-to-entropy ratio and to ob¬ 
tain simple analytic relations among the parameters of 
the model. In this approximation, we consider the co¬ 
moving jV 3 number density to be constant until the time 
^decay = 1 /Tjv 3 , after which it is zero. The N 3 decays 
dump energy into the plasm^f/jand produce scalars that, 
if the plasma temperature at the time of decay is less 
than T c , form mesinos that oscillate and decay on time 
scales short compared to the expansion of the Universe^] 
creating a net baryon asymmetry. We now describe the 
evolution of the Universe in the context of this approxi¬ 
mation. 

At very early times, the Universe’s energy density is 


dominated by radiation, until some time t eq < t, 


decay 


when the energy density in N 3 (matter) is equal to that 
in radiation. For times t eq < t < tdecay? the Universe is 
matter-dominated. After the decays of N% at tdecay, the 


3 Note that solving Eqs. ( |5l[ ) to < |53[ ) exactly, as we will do be¬ 
low, shows that the plasma temperature does not actually in¬ 
crease during N 3 decay—rather, its cooling rate slows BH- For 
this reason we ignore whether N% decays “reheat” the plasma to 
a temperature above T c , a possibility that would naively spoil 
hadronization. 

4 At temperatures below T c , the </>-(/>* annihilation cross section 
becomes large, 0 (Aq CD ) 1421 . However, the rate for this process 
is tiny compared to the oscillation and decay rates since it is 
suppressed by the small ratio of the N 3 to </> decay rates, Tann oc 


Universe is again radiation dominated. The temperature 
of the plasma at the time of matter-radiation equality, 

^■eq? IS 


T — 

- L eq — 


45C(3) 

27 T 4 sr* (Ti) 


rriN 3 , 


(54) 


with < 7 * (7j) the effective number of relativistic degrees 
of freedom at a temperature Ti ~ mjv 3 when the uni¬ 
verse is very hot and radiation dominated. We will use 
9 * (Ti) = 116.25 which takes into account the SM contri¬ 
butions as well as contributions from <j> and IVp 2 - The 
time vs. temperature relationship during radiation dom¬ 
ination fixes the time at T eq , 


^eq — 


45 


Mpi 


2ttV (Teq) T 2 q 


/ 2-7T 2 

3/2 

7 r 2 

V44 

) 

LC(3)J 


gl{Ti) Mp\ 


1/2 

g* 


(Teq) m 


iV, 


(55) 

(56) 


where Mpi = ( 87 rG ) -1 ^ 2 = 2.4 x 10 18 GeV is the Planck 
mass and g * (T eq ) ~ 70, for T eq of order a few GeV. 
After t eq , assuming t /§> t eq the N 3 and radiation energy 
densities evolve according to matter domination, 


_ 4M& ft eq \ 2/3 _ 4Mpj 1 

Pr&d g j , PN 3 - 3 f2 - 


(57) 


Defining £ = p Ns (t = t decay )/p lad {t = t decay ), where 
^decay indicates an infinitesimal time before decay, we find 


5 = 


^decay 

Gq 

45 

27 T 2 


2/3 


C (3) 


4/3 


9* (Teq) 

9t (Ti) 


The baryon-to-entropy ratio is 


Vb = 


n B (t = t 


1/3 


+ i 

decay J 


Srad — t, 


m 


N 3 


\Mp\T N3 


decay1 


(58) 

2/3 

(59) 

(60) 


where t decay denotes an infinitesimal time after decay. 
The net baryon number density is related to the N 3 num¬ 
ber density just before decay, 

n B (* = tdecay) = n N 3 (* = *decay) X ~A e B . (61) 

A and e B are determined solely by the mesino particle 
physics and factor out of the cosmological story. We can 
write B in terms of £, 
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The ratio of the temperatures before and after decay can 
be found by using the conservation of energy, 


Prad ^)decay) Prad ^decay) + P IV3 ^decay 
= Prad (^decay) 6 + 6 > 

which leads to 

3 ^rad (^decay) £ Ae S 

713 4 mjv 3 (l + £) 3/4 2 


(64) 

(65) 


( 66 ) 


We can use the expression for the radiation energy den¬ 
sity in Eq. (57) to find the temperature of the plasma 

just before decay, T ra d (t^ 


T, 


decay — 


— (—) 
Vs \C (3) / 


1/3 r 


15x/3C (3) 1 


decay) ~ +decay, 



g* (Ti ) 

1/4 

f M Pi r A5 

V ' 3 

_ 9* (-^decay ) . 


V ™,n 3 

) 

(67) 

rriN 3 



(68) 

) z ’ 




with g* (Tdecay) = 63.75 the number of relativistic de¬ 
grees of freedom at this temperature. We ignore here 
the decrease the number of degrees of freedom below 
the QCD confinement temperature since it obscures the 
physics and is numerically unimportant at the level of 
accuracy of our estimates. Using this expression, we can 
express the baryon-to-entropy ratio as 


Vb -- 


45 VK (3) 1 


8 tt 4 g * (Tj) (i + £) 3 / 4 


Ae f 


=6.1 x 10” 10 116,25 
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VP* (2i) / \ 1 + £ 

(^)(iF*)- 


3/4 


The factor 


(i + 6 


-3/4 


T r ad (t. 


- V 3 

decay I 


T r ad (t, 


decay1 


(69) 

(70) 

(71) 

(72) 


can be thought of as the entropy dilution of the naive 
i-entropy ratio due to the N :i decays heating the 

For a fixed value of rri m 3 , there is a minimal amount 
of entropy dilution, i.e. a minimal value of £, that can 
be found by setting the temperature of the plasma at 
N 3 decay to be equal to the QCD confinement tempera¬ 
ture. Using this value of the entropy dilution determines 


baryon-t( 

plasmaF] 


5 


Again, we stress that the plasma does not actually heat up in an 
exact treatment of Eqs. (15 1 1 to 1531). See below. 



FIG. 4. Solid lines indicate values of Fjv 3 and mjv 3 that give 
r/ B = 8.6 x 10 -11 for t B = 10 —4 , 10 —5 , and 10 -6 , from bottom 
to top, respectively. The shaded regions lead to a plasma 
temperature at IV 3 decay greater than T c ~ 200 MeV, in 
which case mesinos do not form, or a temperature at decay 
less than 1 MeV, which would spoil BBN. Note that this is 
calculated within the sudden decay approximation. 


a lower bound on the asymmetry e B as a function of mjv 3 
such that a baryon-to-entropy ratio of 9 x 10 -11 is pos¬ 
sible. We plot this lower bound on eg in Fig. [5j as a 
function of the scalar mass, taking the splitting between 
N 3 and (j) to be A m^ 3 = 3 GeV. 

To assess whether we can reasonably achieve a baryon- 
to-entropy ratio in line with the current experimental 
value of 8.6 x 10 -11 , we also show in Fig. [ 5 ] the asymme¬ 
try eg as a function of , taking A= 
1 GeV, y s 1 = y s2 = 10 -2 > y s3 , for a B = 10 -4 , 10” 5 , 
and 10 -6 . We see that there exist ranges of parame¬ 
ter choices that allow for scalars with a mass between 
300 GeV and 10 TeV to have an asymmetry above the 
bound, and therefore to account for the measured value of 
r] B , depending on the choice of Tjv 3 . In principle, masses 
as large as 10 6 ' GeV can be compatible with the ob¬ 
served baryon asymmetry but scales this large would re¬ 
quire a larger e Bl hence a larger amount of fine-tuning. 
In addition, testing this model and embedding it in a 
framework that explains the weak scale, like the super- 
symmetric SM, could become difficult. 

We can check the accuracy of our sudden decay ap¬ 
proximation by solving Eqs. (511 to (53) numerically 
and determining the resulting baryon-to-entropy ratio. 
In Fig. [6j we show the results of a numerical solution 
of the evolution equations with itin 3 = 650 GeV and 
Tn 3 = 10~ 20 GeV. The top panel shows the evolution 
of the plasma temperature in time. As mentioned above, 
the plasma is not “heated up” by the decay of the N 3 s; 
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FIG. 5. The shaded region shows values of tB and m as¬ 
suming = 3 GeV, that result in gg < 8.6 x 10 -11 for 

Tdecay < T c ~ 200 MeV in the sudden decay approxima¬ 
tion. The points above this region can attain the measured 
value of r/B, depending on the value of Vn 3 - Also shown are 
the values of tB when choosing A^jvj = A ( j > N 2 = 1 GeV, 
Vsi = Us 2 = 10~ 2 > y a3 , and a B = 10 -4 , 10 —5 , 10 -6 . 


t [GeV -1 ] 

10 18 10 19 10 2 ° 10 21 



instead, it cools more slowly. The next panel shows the 
ratio of the N 3 energy density to the radiation energy 
density, PAq/pracb as a function of time. This ratio grows 
until ~ 1 /r,v 3 when it decreases to to N :i decays. The 
bottom panel shows the net baryon-to-entropy ratio, gs 
as it develops from zero at early times to a fixed value 
at late times for tB = 10 -4 , 10 -5 , and 10 -6 . We also 
show the value of r]B calculated for these parameters in 
the sudden decay approximation. We note that the sud¬ 
den decay approximation of tjb agrees with the numerical 
solution of Eqs. (51) to (53) to within an ~ O (1) factor 
for a wide range of values of m.N 3 , Tjv 3 , and ts- In our 
numerical solution, we took g * (T) to be fixed at 63.75 
for temperatures less than 4 GeV. Including the change 
of g* at the hadronization temperature would not drasti¬ 
cally affect our conclusions, and would only obscure the 
relevant physics. 


We note here that the estimate of the baryon asymme¬ 
try in this section depended on the assumption that the 
scalars were produced by the late, out-of-equilibrium de¬ 
cay of a particle that was in thermal equilibrium at very 
high temperatures. Alternative cosmological histories are 
certainly possible. For example, the scalars could be pro¬ 
duced by the decay of long-lived moduli [5]. Another pos¬ 
sibility involves considering a nonthermal parent particle 
(possibly N 3 ) that dominates the energy density of the 
universe. Its decays to scalars reheat the universe and 
create a baryon asymmetry before the creation of a ther- 


FIG. 6. Results of a numerical solution of Eqs. ([51| to (531 
for mjv 3 = 650 GeV and = 10 -20 GeV. Top: the temper¬ 


ature of the plasma as a function of time (solid, black). For 
comparison, the evolution of the temperature without N 3 de¬ 
cays, assuming a radiation-dominated Universe is also shown 
(dashed, gray). The dotted blue line shows T c = 200 MeV. 
As mentioned, the plasma does not actually heat up due to N 3 
decays, rather it just cools more slowly. Middle: the evolution 
of the N 3 energy density to that of the plasma. Bottom: the 
growth of r/B for values of the asymmetry tB = 10 -4 (solid, 
black), 10 -5 (dashed, black), 10 -6 (dotted, black) along with 
the value as measured by Planck (solid, red). Also shown are 
the sudden decay approximations to r/B for each of these val¬ 
ues of tg (gray). The decrease in from its maximum to its 
value at late times in each case reflects the entropy dilution 
due to the energy dumped into the plasma. 


mal plasma. Alternative scenarios may, of course, have a 
quite different relationship between ts and ijb from the 
simple one presented here. 


V. CONCLUSIONS 

We have shown that if there is a scalar quark which 
lives long enough to hadronize, CP violation in the os¬ 
cillations of mesinos is possible. If the scalar quark has 
baryon-number-violating decays, and if it can be pro- 
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duced out of thermal equilibrium in the early universe, 
such CP violation could be the origin of the matter- 
antimatter asymmetry of our universe. We presented a 
mechanism for the out-of-equilibrium production of the 
scalar quark via the late decays of a very weakly cou¬ 
pled heavy neutral particle, and computed r]B, the ratio 
of baryon number to entropy, in terms of the relevant 
parameters of the model. Because r)B depends sensi¬ 
tively on the parameters, the viable parameter space of 
the model is rather well constrained. In particular, ob¬ 
taining sufficient CP violation requires the existence of a 
neutral Majorana fermion which is not much lighter than 
the squark. Such a difficult to explain coincidence hints 
at an environmental selection principle, which is attrac¬ 
tive in light of the observation that the value of tjb in 
our universe appears to be optimal for creating stars in 
galaxies which are not so dense as to be unfriendly to 
advanced life [43] . 

In order to experimentally test the model one should 
begin by finding evidence for scalar quarks at the LHC. 
As discussed in Sec. |III A[ existing searches, especially 
in the case where the decays do not produce significant 
missing energy, are not very constraining. In our model 
the squark should decay sometimes into two antiquarks, 
and sometimes into a quark and a neutral fermion, with 
the latter decaying into one soft and two hard quarks. 
The final states are likely to contain at least one third 
generation quark for the reasons discussed in Sec. |IIIB| 
Because successful baryogenesis requires that the neu¬ 
tral fermion not be much lighter than the mesino, the 
first quark is very soft. Thus the scalar quark should 
essentially decay into two hard jets (and possibly two 
very soft quarks). The decays may be prompt or could 
produce displaced vertices. The LHC reach for the dis¬ 
covery of heavy scalar quarks which decay into jets is 
projected to reach 1070 GeV [25] with 1000 fb -1 . In 
Ref. |T2j, Berger, Csaki, Grossman, and Heidenreich pro¬ 
posed looking for same-sign top quarks as evidence for 
mesino oscillations, assuming that the squark predomi¬ 
nantly decays into a top and a strange quark. In our 
model, same sign tops can be produced without mesino 
oscillations, as the neutral fermion which is often pro¬ 
duced in </) decays is equally likely to decay to top and to 
antitop. Searching for such events containing same sign 
tops should increase the reach for searching for scalar 
quarks and significantly constrain the model. Note that 
this signature will be CP-violating, but the small asym¬ 
metry between top and antitop will be too challenging to 
measure at the LHC. 

We also point out that a slight extension to the model 
would allow for asymmetric dark matter (for a review, 
see, e.g., [44]) to be incorporated. If the operator medi¬ 
ating </> decays to baryons, aij<p*diUj, were modified to 
the dimension-5 operator aij<fi*diUj(x/M) where \ is a 
complex scalar and M is the relevant scale, the model 
would possess a Z 2 symmetry under which <j>, and Ni 
are odd. This would render \ stable if it were the lightest 
of these particles. Then the oscillations and decays of the 


mesinos would produce a X~X* asymmetry along with the 
baryon-antibaryon asymmetry. The x mass in this model 
would be fixed to be ~ 5 GeV since the number densities 
of dark and baryonic matter would be equal. The sig¬ 
nature of <j> production at a hadron collider would now 
include missing energy, pushing up the lower bound on 
the allowed mass. 


Appendix A: Form Factor Estimation 

In this appendix, we detail our estimation of the 
mesino-meson transition form factor which enters the 
two-body contribution to T^. We specialize to the 
$ s —>■ r]Ni case that we considered above. 

We begin by expressing the form factor using an opera¬ 
tor product expansion (see, e.g., Ref. [451 for a description 
of the formalism), 

F (d 2 )= [ dxdy(^{p')\^ a s^\0)C l ^ s (0\(t^s l s \^ a p (p)). 

(Al) 

p is the mesino momentum and p' = p + q is the me¬ 
son momentum—a transition involving an on-shell N\ 
has q 2 = m^. is the short-distance Wilson co¬ 

efficient calculated perturbatively. Greek subscripts are 
Lorentz indices and Latin superscripts are color indices. 
We use the shorthand f dx = J 0 ' dx\ f ( ' dx 2 S(l — x\ — X 2 ) 
and similarly for dy. X\ is the momentum fraction of the 
squark in the mesino and X 2 is that of the light quark. 
1 / 1,2 are the momentum fractions of the quark and anti¬ 
quark in the meson. We expand the bilinear operators 
in terms of operators of definite quantum numbers using 
the meson and mesino momentum distributions]^] 

{v{p')Ksj\ 0 } = i / 7 5 ) QT ^ ( 2 / 1 , 2 / 2 ) + 

(A2) 

, S jl 

(0|<^s 5 |4> s/3 (p)) = —8p S 4>^. (x 1 ,x 2 ) + ■■■ , (A3) 

where the distribution amplitudes </>,, and <j)$ s are nor¬ 
malized so that 

/**,<*,») = /„ J dx<j)# s (x i,x 2 ) = U a , (A4) 

where f v and /$ 3 are the meson and mesino decay con¬ 
stants. 

At leading order there are two diagrams that con¬ 
tribute to the Wilson coefficient, shown in Fig. [TJ which 


We ignore the uu and dd quark components of 77 , an approxima¬ 
tion that is valid given the overall uncertainty of the calculation. 
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-O 


g' 



o 


momentum fraction m g /m$ s where m q ~ Aqcd is a con¬ 
stituent quark mass. 

Integrating over the mesino momentum distribution 
takes 


FIG. 7. Leading diagrams responsible for the $ s —> r/Ni 
transition. Dashed lines are <j> lines and solid lines are s or 
s lines. Curly lines are gluons. Shaded dots represent the 
tfi-s-Ni coupling. 


give 


C i ^s = (*9sf(-^)h P l s (tT(n kl 


k 2 


A? 


Pr - 


IPJ'yS ' 

i(x 1 p + A 2 ) p 
-ml 


Pt 


(A5) 


a/3 


Ax = p’ — x 2 p and A 2 = p — y 2 p' are the momenta of 
the virtual quark and squark in the two diagrams and 
k is the gluon momentum. Using these, the form factor 
becomes 


V • Ax 

A? 

V • (xxp+ Ax) 
A 1 -m? 


-2 + 0 


Aqcd 


2 2 
to| - p 


m% s — — y 2 (m% s - q 2 ) 


■O 


Aqcd 

m$ s 


and 

fc 2 ~ -AQCDTO$ s y2 ( 1 - 


TO 3 


■ 0 


A 2 




(All) 

(A12) 


(A13) 


Because of the factor of 1/fc 2 in the Wilson coefficient, the 
form factor is enhanced in the region where q 2 approaches 
to| giving the form factor a 1/ (rn% — q 2 ) scaling. Af¬ 
ter integrating over the meson momentum distribution, 
(A121 gives a logarithmic correction to the scaling which 


we can ignore. We therefore take C ~ —2, finding 


n<f) 


= 7r a s 


tr (t a t a ) 
N* 

1 *c 


(A6) 


x 


J dx dy (j> M (xi,x 2 )4>* P (?/i, y 2 ) 


C_ 


F(d 2 ) 


3 (A 2 -1) /„/*. 1 

—-- -7T(y ----- 

N 2 s A qcd to$ s 1 - q 2 /m% 
87ra s f v f$ s _ 1 

3 Aq C dto$ s 1 - q 2 /m% s ' 


(A14) 

(A15) 


with 


C 1 = ^2 Tr £/ 7 V4 i P Llp ) (A7) 

- a2 * 2 Tr |/ 7 5 P l (x\j/+ ^i)) 

^2 m <j> 

= tUAi + y-(*ip+A,) 

A 2 A 2 -to 2 

We will use simple forms for the distribution functions, 

07? ( 2 / 1 , 2 / 2 ) = 6/^ 2 / 13 / 2 , (A9) 

( 31 , 3 : 2 ) = UJ (x 2 - Aqce> N ) . (A10) 

V ) 

where Aqcd ~ 300 MeV is the QCD scale. The me¬ 
son distribution function is motivated by the requirement 
that the quark and antiquark momentum distributions be 
symmetric—it can also be shown to be the asymptotic 
meson distribution amplitude. The mesino distribution 
function encodes the fact that the light quark carries a 


Using a s ~ 0.3, appropriate for the mass splitting be¬ 
tween and Ni of ~ 1 GeV that we have in mind, 
Aqcd — 300 MeV, a value of the rj decay constant de¬ 
rived from the 77 —> 77 rate (f v = 130 MeV), and the 
mesino decay constant from Eq. (21 1 , this gives, for the 
production of a physical N \, 


F(m 2 Nl ) ~ 1.2 x 10 


-2 


1 650 GeV ( 1 GeV 


m$> s 


Am 


<E>iVi 


(Aie) 


Being an estimate of an exclusive process in QCD, this 
calculation has considerable uncertainty, and represents 
the largest contribution to the overall uncertainty in our 
estimate of the baryon-to-entropy ratio. However, it cap¬ 
tures the essential physics: this rate is enhanced when the 
singlet mass is close to the mass of the mesino. This is 
because the quark produced in the (f> -+ Nis splitting 
has a smaller momentum in this regime, allowing it to 
more easily bind with the spectator quark to form a me¬ 
son. Similar behavior is observed in the B meson system, 
where the semileptonic transition form factor in B —> ntv 
decays grows as the tv mass approaches ms- 

Note that the overall normalization of the form factor, 
for a fixed |Ti 2 1, is degenerate with the value of \y s i |. 
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